Conformal, Subconformal and Spectral Universality in Incommensurate Spin Chains 
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A renormalization scheme is developed to study an 
anisotropic quantum XY spin chain in a quasiperiodic trans- 
verse field. The critical phase of the quasi-particle excitations 
of the model with fractal wave functions exists in a finite pa- 
rameter interval and is sandwiched between the extended and 
localized phases. The scaling properties of the critical phase 
fall into four distinct universality classes referred as spectral, 
subconformal, conformal and Harper. The spectral and con- 
formal classes respectively describe the onsets of extended to 
critical and critical to localized transitions while the subcon- 
formal class describes the part of the phase diagram sand- 
wiched between the conformal and spectral transitions. The 
Harper universality class describes the isotropic limit of the 
XY model. A decimation scheme is developed to compute the 
infinite sets of universal scaling ratios characterizing the wave 
functions in the four universality classes. The renormalization 
flow equations exhibit a limit cycle at the band center and at 
the band edges providing a new method for determining these 
energies with extremely high precision. 
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I. INTRODUCTION 

In recent years, renormalization group (RG) techniques 
have been the dominant theme in describing the scal- 
ing properties at the onset of phase transitions. These 
studies include the thermodynamical phase transitions, 
driven by thermal fluctuations, as the temperature is var- 
ied as well as the transitions which occur when a control 
parameter is varied. The later include the transitions to 
chaos in nonlinear dynamical systems and the phase tran- 
sitions in quantum systems at zero temperature such as 
the metal- insulator transitions, driven by quantum fluc- 
tuations. 

One-dimensional quantum spin chains provide very 
simple and interesting models that exhibit phase tran- 
sitions at zero temperature as a function of anisotropy or 
magnetic field. ||^ In some cases, the T = phase tran- 
sitions of the quantum spin chain map onto the phase 
transitions of a 2D classical model in which the vari- 



able parameter is the temperature. These systems are 
also of experimental interest as various physical realiza- 
tions of quantum spin chains have been found in quasi- ID 
magnetic systems. [D The Ising and the anisotropic XY 
model in a constant transverse field provide the simplest 
class of quantum models [|| (that can be solved exactly) 
that exhibit the second order magnetic phase transition 
at a critical value of the magnetic field. The anisotropic 
model belongs to the universality class of the Ising model, 
i.e. the critical exponents at the onset of transition are 
the same as those of the Ising model ||]. Therefore, such 
anisotropic models with a uniaxial anisotropy are usually 
referred as Ising like. 

Motivated by the discovery of a quasi-crystal and 
the fabrication of magnetic superlattices, quantum spin 
chains in a quasiperiodic (QP) transverse field have been 
investigated extensively. Initial studies were confined to 
the case where the transverse magnetic field was inho- 
mogenious with binary Fibonacci distribution. Q| The 
quasiparticle excitation spectrum turned out to be a frac- 
tal for all parameter values. These studies were fol- 
lowed by other ones in which the magnetic field was os- 
cillatory with the period incommensurate with respect 
to the period of the spin chain Detailed numeri- 
cal work showed that the quasiparticle excitations of the 
anisotropic quantum spin chain could be extended (E), 
localized (L) or critical (C) resulting in E-C and C-L 
transitions as the strength of the magnetic field was var- 
ied. These three different phases were respectively char- 
acterized by continuous, point like and singular contin- 
uous spectrum. It should be noted that in the isotropic 
limit the eigenvalue equation for the quasiparticle exci- 
tations of the model reduce to the the famous Harper 
equation exhibiting both extended and localized phases 
and hence the metal-insulator transition. |^ The exis- 
tence of a self-similar butterfly spectrum at the onset of 
transition aroused a great deal of interest in this model. 
The Harper equation has also been of interest due to 
the fact that it describes the two-dimensional electron gas 
in a magnetic field. Q Furthermore, this model has at- 
tracted the attention of mathematicians due to the small 
denominator problems. |^] 
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The anisotropic XY model describes a perturbed 
Harper model and therefore provides a new class of 
quasiperiodic models which splits the single E-L transi- 
tion of the Harper model into two transitions, namely 
E-C and C-L transitions. The existence of a fat C 
phase sandwiched between E and L phases provides a 
new paradigm for the metal-insulator transition in one- 
dimensional incommensurate systems. Analogous to the 
Harper model, in the anisotropic model both the E-C 
and C-L transitions are global in the sense that all quan- 
tum states are E, C or L simultaneously. Fig. 1 shows 
the spectral phase diagram of the model. The line in 
the phase diagram corresponding to the onset of local- 
ization also corresponds to the onset of magnetic phase 
transition from ordered to disordered phase providing an 
interesting example of a spectral and magnetic interplay. 
This line will be referred as the conformal line as along 
this line, the gap in the excitation spectrum vanishes and 
the model is conformally invariant. The E-C transi- 
tion in the quasiparticle wave functions does not affect 
the magnetic properties and hence the corresponding line 
will be referred as the spectral line. The fat C region 
sandwiched between the the spectral and conformal line 
will be called subconformal. The conformal and spectral 
lines intersect at the Harper limit. 

In our previous short paper [Tl[ | , we described a new 
decimation scheme to study nearest-neighbor (n.n.) tight 
binding models (TBM) and applied it to QP Ising chains 
exhibiting the C-L transition. The proposed scheme con- 
firmed the existence of a fat C phase in the Ising limit. 
It also demonstrated the fact that the fat C phase was 
described by two distinct fixed points which respectively 
described the C-L transition at the conformally invariant 
point and the subconformal C phase. Since the deci- 
mation scheme was valid only for n.n. TBM's, it de- 
scribed only the Ising limit of the general anisotropic 
model containing also a next-nearest-neighbor (n.n.n.) 
term. Therefore, these studies did not provide a full char- 
acterization of the C phase. In particular, the universal- 
ity class for the E-C transition and for the conformal line 
(see Fig. 1) remained an open question. 

In this paper we generalize our method to study the 
anisotropic model which includes Ising as a special case. 
The scalar decimation equations of our previous studies 
are replaced by vector equations. We show that the C 
phase is described by four distinct universality classes. 

In section 2, we describe the model and the corre- 
sponding tight binding form. In section 3, we discuss the 
decimation scheme based on the Fibonacci lattice sites. 
Sections 4 and 5 respectively deal with the qualitative 
as well as quantitative characterization of universal fea- 
tures of the phase diagram. In section 6, we show how 
the decimation scheme can be generalized for an arbitrary 
path in the Farey tree. In section 7, we discuss the RG 
equations. Section 8 summarizes the various universality 
classes, and we end with some conclusions in Section 9. 



II. SPIN HAMILTONIAN AND THE TIGHT 
BINDING EQUATION 

The anisotropic XY model in a transverse field is given 

by 

^ = - Et^^-^-^+i + Jy'^^'^hi) + ''^'^n (1) 

i 

Here a are the Pauli matrices describing the spin 1/2. 
Jx and Jy describe ferromagnetic exchange interactions. 
In our discussion, unless otherwise stated, hi is an in- 
homogcnious modulating field of periodicity a which is 
incommensurate with the periodicity of the spin chain: 

= Xcos{2Tr{ia + (f>)) (2) 

We choose the parameter a to be equal to the inverse 
golden mean a = (-\/5 — l)/2. (/)isa constant phase 
factor. 

Using the Jordan- Wigner transformation, the model 
can be transformed into that of a free spinless fermion 
Hamiltonian 

+ hic\ci h.c], (3) 

We see that the lack of the 0{2) symmetry in the spin 
Hamiltonian results in a lack of the U{1) symmetry 
(e*^Ci — UciW) in the corresponding fermion Hamilto- 
nian. The system can be diagonalized using the methods 
described by Lieb et al. ||^. The resulting coupled eigen- 
value equations for the quasi-particle excitations are |^ 

Jyipi+i + Jxipi-i + 'ihiipt = Erii (4) 
JxVi+i + JyVt-i + '^h^m = Eil>i (5) 

The above two coupled equations with n.n. terms can be 
combined into a single equation with n.n.n. interaction: 

+ 2(Jy/ii_i -I- Jxhi)^.,_i + 2{Jyhi + Jxhi+i)ipi+i 
+ {J^ + + 4hf)i;„ (6) 

Alternatively, the coupled equations can be rewritten as 
a single vector equation involving only a n.n. interaction: 

J{y, x)^i+i + J(x, + 2/1,*, = E'fi, (7) 

where J{x, y) is a 2 x 2 matrix whose nondiagonal entries 
are zero while its diagonal elements are Jy and Jx (see 
Eqs. 4 and 5) and ^' is a two-dimensional vector. 

In the isotropic limit, the coupled equations become 
degenerate and the two components of ^> become identi- 
cal. The resulting equation is the famous Harper equa- 
tion 1^ which is a paradigm for quasiperiodic models ex- 
hibiting the metal-insulator transition. The anisotropic 
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model describing a perturbed Harper equation can also 
be written as 



(8) 



where Hh denotes the Harper Hamiltonian consisting of 
the first and the third term in Eq. (3) and Hg denotes 
the part of the Hamiltonian which exists only at a non- 
zero value of the anisotropy g defined as Jx = 1 — g and 
Jy = 1 + g. We see that the anisotropy perturbs the 
Harper squared model. 

The fact that the anisotropy doubles the size of the 
Hilbert space can be seen also in the semiclassical limit. 
The anisotropic model in the semiclassical limit describes 
the motion of a particle with the Hamiltonian 

1 i 
H = -{Jx + Jy)cos(p) + -{Jx - Jy)sin{p) + \cos{x), 

(9) 

where P = 7^ is the momentum operator. We see that 
the anisotropy perturbs the system along a complex di- 
rection. This explains why the Hilbert space has doubled 
its dimension in the anisotropic case. 

In the Ising limit Jy = 0, the matrices J(a;, y) of the 
vector equation are noninvertible. In this case, we will al- 
ways work with the TBM which contains only n.n. terms: 

2Jxh,+,ij,+i + 2Jxh,ij,_i + {Jl + 4/i2)Vj, = E''^,. (10) 

An interesting limit of the Ising model is the case Jx — *■ 
00 and A — > ||l^ where the TBM contains no diagonal 
disorder: 



2hi+iipi+i + 2hi%pi. 



(11) 



It is easy to show that the model is always self-dual and 
has been numerically shown to have only critical states. 
Therefore we refer to it as the critical model. It should 
be noted that this TBM also describes the isotropic XY 
model with oscillatory exchange interaction Jx(i) — hi 
and a constant field. |13| This should be contrasted with 
the Harper model which describes the isotropic XY model 
with constant exchange and oscillotory field. 

The numerical phase diagram 1^ for the anisotropic 
model describing a perturbed Harper model is shown in 
Fig. 1. Whenever A and one of the exchange inter- 
actions become equal either the E-C or C-L transition 
takes place. The E-C transition corresponds always to a 
smaller absolute value of A, i.e. to the exchange interac- 
tion whose absolute value is smaller. Whenever Jx and 
Jy differ a fat C phase is observed in the phase diagram. 
Therefore, the fattening of the critical point is due to the 
breaking of the U{1) symmetry which is a consequence of 
the broken 0{2) symmetry in spin space. An interesting 
consequence of the spectral and magnetic interplay is the 
fact that the onset of the C-L transition is coincident with 



the onset of the magnetic transition to the long range or- 
der (LRO) where two-point long range spin-spin correla- 
tions vanish. These numerical observations regarding the 
onset of E-C and C-L transitions and the fact that C-L 
transition coincides with magnetic transition still remain 
unexplained. 

The energy spectrum of the anisotropic model is char- 
acterized by a gap in the spectrum which vanishes along 
the conformal line. Therefore, unlike the isotropic limit, 
_B = is not an eigenvalue of the system in most part 
of the phase diagram. As seen below, this introduces an- 
other degree of complexity in analysing the scaling prop- 
erties of the model. 



III. FIBONACCI DECIMATION 
III. A Formulation 

One of the challenging questions posed by the fat C 
phase is how the scaling properties vary within the phase. 
A priori there are two main possibilities: either the scal- 
ing can be described by only few different universal- 
ity classes or the scaling properties vary continuously 
through the phase. The former case would correspond 
to a few distinct RG fixed points whereas the latter case 
could be described by a line of RG fixed points in function 
space. 

In recent years, various RG schemes have been pro- 
posed to study the scaling properties of the critical point 
of the Harper equation, jit), Our methol- 

ogy is somewhat similar to that of Ostlund et al. ||l^ as 



it describes the scaling properties of the wave functions 
for specific values of the energy. Furthermore, in analogy 
with their scheme, our decimation scheme has been ap- 
plied to study the breakup of KAM tori in the circle and 
standard maps. |jl^, However, the method of Ostlund 
and Pandit jlj], which was based on transfer matrices, 
had limited success due to the fact that at the localization 
threshold infinite products of transfer matrices diverged 
for almost all values of the phase 4>. We propose a new 
decimation scheme where instead of multiplying transfer 
matrices, the TBM itself is decimated. The main advan- 
tage is the reduction in the number of functions needed 
to carry out the renormalization. The cost we have to 
pay is that our recursion relations will be slighly more 
complicated. However, it turns out that with fewer func- 
tions we are able to eliminate directions which lead to 
divergences. This not only helps in approaching the RG 
problem but also provides practical means of calculating 
various essential quantities like the localization threshold 
and the minimum {Emin) and maximum (Emax) energy 
eigenvalues in general TBM's. The knowledge of ma- 
chine precision Emm and E^ax comes from the study of 
the states with Bloch index k = and k — 1/2 in addi- 
tion to fc = 1/4 which respectively correspond to Emim 
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Emax and E = quantum states. The previous stud- 
ies of Ostlund and Pandit and those of Dominguez, 



Wiecko and Jose |17| were restricted to the Bloch index 
/c = 1/4 corresponding to E = part of the spectrum. 

We will first describe the decimation scheme based on 
the Fibonacci lattice sites. This procedure describes only 
three critical scahng ratios out of infinity that are re- 
quired to fuUy characterize the wave function. However, 
the Fibonacci decimation scheme is quite sufficient to 
confirm the fat C phase and to calculate critical param- 
eters and energies. A decimation scheme based on an 
arbitrary path in the Farey tree (see sections 6 and 7) 
will determine the infinite set of universal scaling ratios 
needed to describe the universality completely. As dis- 
cussed later in the paper, the dominant peaks of the wave 
functions are described by decimation. 

In this scheme, all sites except those labelled by the 
Fibonacci numbers Fn {n = 0,1,2,...) are decimated. 
This results in a TBM connecting a function "if at two 
neighboring Fibonacci sites: 



+ Fn+l) = c„(i)*(i + Fn) + d„(i)^'(i) 



(12) 



In the Harper and the Ising limits, the above equation 
is a scalar equation with scalar "decimation" functions 
c„(i) and (i„(i). For the general anisotropic case, ^'(i) 
is taken to be the vector (jpi,r]i) of the coupled eigen- 
value equation and the multiplying factors c„ and dn are 
2x2 "decimation" matrices. In this case, the name "dec- 
imation function" refers to the entries of these matrices, 
which are functions of the lattice site i. The index n la- 
bels the " level" of the decimation. As the fractional part 
of Fna decays to zero as n ^ cxd, we expect the scal- 
ing of the wave function for consecutive Fibonacci sites 
to show some regularities which should also come up in 
the decimation functions. Using the defining property of 
the Fibonacci numbers, Fn+i = Fn + Fn-i^ the following 
recursion relations are obtained for c„ and dn fl^ : 

Cn+lii) = c„(i -I- Fn)Cn-lii + Fn) - ^ (i)d„+i (i) (13) 
dn+lii) = ~dn{i)[dn{i + Fn) + 

Cn{i + Fn)dn-l{l + Fn)]c-\j) . (14) 

As seen from the above equations, the decimated ma- 
trices describe a flow in Fibonacci space provided their 
inverse exist. 

In order to iterate these equations, we need the initial 
conditions for the matrices at the levels n = 1 and n = 2, 
say. For the standard Fibonacci numbers with Fi ^ F2 = 
1, it is possible to take ci to be the unity matrix and let 
di vanish. The non-trivial initial conditions at n = 2 are 
obtained directly from the defining equations (4-5) for 
the TBM. The decimation could be carried out also with 
other choices for Fi and F2 but this would require more 
elaborate initial conditions. 

Comparing our decimation scheme with that of 
Ostlund et al. [p3 , it should be noticed that our approach 



based on the special TBM form requires half the num- 
ber of functions needed for the method based on transfer 
matrices. As discussed below, our new approach elimi- 
nates most of the divergences encountered in the previous 
approach. 

III.B Results and applications 

The recursion relations for the decimation functions 
were iterated numerically for extremely large size systems 
(upto sizes 500,000) required to study lattice of F24 spins. 
The renormalized 2x2 decimation matrices were then 
diagonalized. 

Numerical iterations of the scalar decimation equa- 
tions for the Harper and the Ising models result in well- 
defined asymptotic solutions for the functions c„ and d„ 
as n — > cx) in both E and C phases. For the vector deci- 
mation of the general anisotropic model, asymptotic dec- 
imation matrices and their inverse were found to exist in 
the E phase as well as along the E — C transition line 
in the phase diagram. However, along the conformal line 
and within the fat C phase for the general anisotropic 
case (except the Ising limit), the c— or d- matrices were 
non-invertible for some lattice sites and hence our deci- 
mation equations could not be used there. We used al- 
ternative means to study this part of the phase diagram 
as described in the next section. 

The asymptotic behaviour of the decimation functions 
as n — > 00 is different in the E,C, and L phases and also 
depends sensitively on the energy E and the phase (j). 
In the L phase, c„ always diverges and dn tends to zero 
as one would expect by looking at Eq. (12). In the E 
and C phases, the decimation functions remain bounded 
for all n (except for some cases in the E phase, see Ta- 
ble I). Furthermore, at the band edges corresponding to 
the minimum and maximum eigenvalue Emm and Emax 
(and also at the band center for Harper), the decimation 
functions converge on well-defined limit cycles provided 
the phase factor is tuned to some special values. For 
arbitrary (j) and also for other eigenvalues, the decimation 
functions oscillate in a rather irregular way with increas- 
ing n converging perhaps on a strange attractor. 

Table 1 shows the trivial limit cycles for the E phase. 
The asymptotic solutions are independent of the chosen 
value of (j) and hence do not depend upon the site index 
i. This provides a rather interesting example of the di- 
mensional reduction for a system with infinite degrees of 
freedom. It should be noted that the limit cycle is almost 
independent of the anisotropy g. 

For the Harper model and also along the conformal line 
in the phase diagram the existence of a six-cycle can be 
seen for the zero-energy state. For rest of the phase dia- 
gram, the model has no zero-energy eigenstate and hence 
the determination of a limit cycle requires eigenencrgy to 
an extremely high precision. The values of the other pa- 
rameters determines how sensitive the solution is to an 
error in the energy. Typically, in order to see the limit 
cycle with three digit precision, energies are required to 
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have machine precision. In the worst case of a delayed 
crossover, where the parameter values are very close but 
not exactly on a transition line, it may happen that ma- 
chine precision energy may give an asymptotic limit cycle 
with single digit precision. Obtaining machine precision 
energies is almost an impossible task even for the tridi- 
agonal matrices. (See Table III) Table II illustrates how 
the fixed point found in the E phase can be used to nu- 
merically determine E^ax for an arbitrary value of A. A 
good initial estimate of E„iax can be obtained via a di- 
rect diagonalization of the TBM with periodic boundary 
conditions. This initial estimate can then be improved 
by imposing the fixed point property. The resulting en- 
ergy eigenvalues are much more accurate than the ones 
obtained by the diagonalization. 

Unlike the trivial limit cycle of the E phase, the C 
phase is characterized by non-trivial asymptotic 6-cycles 
observed at the band edges for carefully chosen values of 
4>. Furthermore, the decimation functions depend upon 
i (see Figs. 5 and 6). Our numerical results show that 
(f) — 1/4, for which the QP potential has the symmetry 
h-i = —hi, always gives a desired six-cycle throughout 
the fat C phase. The value = leads also to a reflection 
symmetry (due to various symmetries in the anisotropic 
model the phase can be defined mod 1/2). The phase 
<f) = 1/4 is however special because for the part of the 
phase diagram where zero energy state is an eigenstate 
(Harper model and the anisotropic model along the con- 
formal line), this value of (j) causes the main peak of the 
wave function to lie at i = and the resulting wave func- 
tion is bounded. (See Section IV) 

In the cases where £' = is not an eigenstate, 0=1/4 
was found to shift the main peak in a very orderly fash- 
ion (related to number theory) for finite lattices. For 
the Fibonacci lattices of size F^, the main peak is lo- 
cated at pk, where pk satisfies the recursion Pk+i — Pu = 
^{Pk - Pk-i) + {Pk-i ~ Pk-2)- Therefore, the shifts in 
the peak positions pk — Pk-i give the rational approx- 
imants of a^. For the odd Fibonacci lattice sizes the 
shifts Pk — Pk-i, are given by the even Fibonacci se- 
quence 8,34,144.... while that for the lattices of size 
8,34,144..., the shifts are 21,89,377, .... The shift in the 
main peak between two neighboring odd-even cases is 
4, 17, 72, 305, .... This systematic shift in the main peak 
implies that the continuously varying critical <j} for locat- 
ing the main peak at i = can be found upto arbitrary 
precision by the formula = limfc^oo < cpfc + 1/4 >, 
where < > denotes the fractional part. For the fat C 
phase, this phase factor is found to vary continuously as 
Pq,Pi,P2 vary inside the fat C phase. The knowlege of (f> 
is essential in confirming the universality as discussed in 
the next two sections. 

Detailed numerical iterations of the decimation equa- 
tions reveal that the fat C phase of the model is described 
by four distinct limit cycles: 

(1) Spectral limit cycle along the E-C transition line. 



(2) Conformal limit cycle along the conformal line (C-L 
line). 

(3) Subconformal limit cycle, for the fat C regime 
bounded by E-C and C-L transitions. 

(4) Harper limit cycle corresponding to the isotropic limit 
of the model. 

As stated earlier, all the four limit cycles were found for 
Emin as well as for E^ax (and also for E=0 for Harper). 
It is interesting that the spectral limit cycle does not 
depend upon the actual value of the Emin and Emax 
which varies continuously along the E-C transition line. 
For the subconformal limit cycle, both the energy and 
the strength of the field A were irrelevant parameters. 
Therefore, the subconformal limit cycle is identical to 
the limit cycle observed for the critical model (Eq. 11). 

Table HI compares the eigenenergies obtained by diag- 
onalizing tri-diagonal matrices of varying sizes with those 
obtained by the Newton method based on the existence 
of the six-cycle. This table clearly shows the superiority 
of the decimation method compared to the exact diago- 
nalization even for tri-diagonal matrices. As stated ear- 
lier, a very precise value of Emin is essential to confirm 
the C phase at subconformal points. Near the confor- 
mally invariant point we observe that Emin behaves like 
Jx^^'cil - Xy with c « 0.3708641926 and z « 1.38897. 
This gives a rather good estimate for Emin through the 
whole fat C phase. Close to the limit A = 0, Emin be- 
haves like Jx + E^^^X/2 + where £^^i„ is the band 
edge of the critical model. This fact together with Eq. 
(11) implies a continuous connection between the band 
edge of the critical model and Emin of the Ising model. 

In order to confirm the universality of the observed 
limit cycles, the decimation procedure was repeated for 
the multiharmonic QP potentials 

h{i) = ■.{cos{2'K{ia + (f)) + acos(67r(io' + 4>))) 

Vl + 

(15) 

The observed limit cycles for this generic potential were 
identical to those observed for a single harmonic field in 
all the universality classes discussed above. Therefore, 
we conclude that the fat C phase persists for generic po- 
tentials. 

Numerics to demonstrate the limit cycle in the deci- 
mation functions for multi-harmonic field is complicated 
by the fact that the phase boundaries have to be known 
to a very high precision. We were again assisted by the 
existence of limit cycles which provided a very efficient 
Newton method to locate the phase boundaries. In Table 
IV, we display the flow of the decimation functions to an 
asymptotic 6-cycle for the Ising model at the onset of the 
C-L transition where the model is conformally invariant. 
The results of this table are used in determing the C-L 
phase boundary as shown in table V. 
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The region of the phase diagram where the vector dec- 
imation cannot be appUed due to noninvertibilty of the 
decimation matrices, will be studied in the later part of 
the next section. Based on those studies, we conjecture 
that these regimes do not introduce any new limit cycles. 



IV. SYMMETRIC AND ASYMMETRIC WAVE 
FUNCTIONS 

The existence of a p-cycle for the decimation functions 
often implies that the self-similarity in the wave function 
is described by the equation, 



«'(i)«^'([crPi + l/2]) 



(16) 



([ ] denotes the integer part). This equality with p ~ 6 
(or p = 3 if the absolute values are considered) was origi- 
nally found for the Harper model at its critical point. Our 
studies generalizes the validity of this result for the whole 
class of TBM's discussed here. In case where both the 
decimation functions and the wave function are bounded, 
a limit cycle for the decimations functions necessiates the 
above-form self-similarity for the wave function and vice 
verse. 
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For a given eigenstate with known eigenvalue E, the 
TBM provides an iterative scheme to obtain the wave 
function at all sites in terms of ^'(0) and ^'(1). We choose 
"00 — 1 and adjust the phase factor so that the maxi- 
mum of the absolute value of the wave function -0 is at 
1 = 0. This method always gives a bounded solution for 
the wave function. The unknowns ryo, V'l, and rji are 
then determined requiring that the above equation be- 
comes exact as i tends to infinity. Eq. (12) implies 



where 



(17) 



(18) 



Cn,p = C„+p_i(0)c„+p_2(0)...C„(0) 

p-1 

Dn,p ^^Cn+i+l,p-t-ldn+i{0) {Cn+pfi = Unity) (19) 
1=0 

If there is an asymptotic p-cycle for the decimation func- 
tions, Gn,p and Dn^p approach a fixed point in the limit 
n ^ oo. Assuming that 4'(F„+p) = ^'(-F„) = C*(0) we 
can write Eq. (17) as an eigenvalue equation and solve 
for the eigenvalue C and the eigenvector *(0) (i.e. 770). 
An estimate for \I'(1) is then obtained by applying the 
recursion relation (12) backwards n + p—2 times. This 
section describes wave functions obtained in this way. 

Figs. 2 and 3 show the wave functions corresponding 
to different parts of the phase diagram. It is interesting 
to note that for Harper, the wave function is symmetri- 
cal about i = Q wherease for the Ising case at the con- 
formally invariant point, it is completely asymmetrical. 



The symmetry in Harper is obvious from the TBM. The 
asymmetry of the E = Q eigenstate in the Ising case is a 
consequence of the fact that in this limit the wave func- 
tion is given by the simple recursion J^tpi-i + 2hi'ipi = 0. 
Because ho — with = 1/4, we obtain ipi = for 
all i < 0. A more physical picture of the wave function 
could be obtained by choosing another phase factor for 
which hi would be nonvanishing for all finite i. However, 
even in this case there would always exist an i so that 
hi K and the wave function would essentially vanish on 
a finite left-hand-side neighborhood of that lattice site. 
The asymmetry can be understood also from the semi- 
classical analysis of the TBM. jl^ In the Ising case, the 
eigenvalue equation in the semiclassical limit is (from Eq. 
9) 

exp{—ip)ijj{x) + \cos{x)ijj{x) — Er/{x) (20) 
exp{ip)ri{x) + \cos{x)ri{x) — Eip{x) (21) 

For E — 0, if cos{xo) — for some ccq, then 
exp{—ip)'ip{xo) = exp{ip)'r]{xo) = 0. Since p is the gener- 
ator of space translation, the above equation implies that 
'ijj{x) and ri{x) respectively vanish to the left and right of 

Xq. 

The figures suggest that the phase diagram for arbi- 
trary g can be classified as either symmetric or asym- 
metric. Along the conformal line in the phase diagram, 
the wave function is always asymmetric about the cen- 
ter. However, in the subconformal region of the C phase 
the system slowly recovers its symmetry asymptotically. 
For the critical model, the wave function fully restores its 
symmetry and in fact qualitatively resembles the Harper 
case. For the critical model, the phase = (3 — 2<t)/4 
setting the main peak at i = is such that the poten- 
tial has the symmetry h^i = —hi^i. This implies the 
reflection symmetry for the wave function. Furthermore, 
the E-C transition is also characterized by an asymptot- 
ically symmetric wave function. Fig. 4 shows the Hull 
function, obtained by plotting the wave function \ipi\ as 
a function of < icr >, in the four universality classes. 

In the region where the vector decimation could not be 
implemented, we have to use other methods. Along the 
conformal line for £' = state, the coupled equations (4- 
5) decouple. It turns out that either the forward or back- 
ward iteration of Eq. (4) is extremely stable beginning 
with an arbitrary initial condition. In other words, we 
can give 4^±n an arbitrary initial value and iterate in the 
stable direction to obtain ipi for i = ...—A/, — Af+1, M. 
These '0-values are very accurate if N is much larger 
than M. In the subconformal region of the anisotropic 
model, the wave function is obtained by the diagonaliza- 
tion method and thus it is there not so accurate as in 
other parts of the phase diagram. 

In the next section, we provide a quantitative charac- 
terization of the universality classes. 
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V. INFINITE SET OF UNIVERSAL SCALING 
RATIOS 

Ostlund and Pandit jl^ described the universality of 
the Harper critical point in terms of the scaling ratio 



c 



lim |^(F3„)/V'(0)|, 

n— >oo 



(22) 



This scaling ratio along with Eq. (16) clearly describes 
the wave function in the C phase by implying that it is 
neither extended nor localized. Eq. (22) describes the 
decay of wave functions at asymptotic Fibonacci sites 
with respect to the central peak. The existence of a limit 
cycle with period 3 in the decimation equation (consid- 
ering the absolute values) gives three different values of 
(. Two of them are listed in Table VI. The interesting 
question is what happens to the similar scaling ratios at 
sites other than the Fibonacci's. Does the wave function 
repeat itself in the same way as it does for the Fibonacci 
sites? Our detailed numerical studies suggest that it is 
true for all sites that a critical wave function repeats it- 
self: i.e for every given site, there exists a whole sequence 
of sites Qk where the wave function approaches the same 
amplitude which is a universal fraction of the main peak. 
This sequence is given by the recursion relation 



Qk+i = 4:Qk + Qk-i- 



(23) 



The recursion is such that Qk/Qk+i ^ cr^ a.s k ~^ oo. 
This implies that the scale invariance of the wave function 
is described in terms of an infinite set of scaling ratios C, 



C(Qo,Qi) 



lim \tP{Qk)/m\, 

k—>-oo 



(24) 



where Qk is obtained by the above recursion relation. 
The set obtained by varying the two integers Qo,Qi is 
complete and hence specifies the location and height of 
all the peaks of the wave functions. Stated differently, 
this result implies that the wave function at every Qk 
site repeats itself at sites given by the above recursion 
relation. It is easy to see that this generalization includes 
the previous results [Q related to the Fibonacci sites as 
a special case. (See Table VI) 

We next investigate the sequences that result in the 
dominant peaks in the wave function. Our detailed stud- 
ies show that out of infinity of the scaling ratios, the C's 
for the dominant peaks are obtained from those values of 
Qo and Qi that are of the order of 4. 

Furthermore, it turns out that the dominant peaks can 
be labeled in two ways. This also helps in regrouping the 
infinite set of s. First, with every sequence determined 
by Qo and Qi, we associate a set of n "harmonic" se- 
quences {nQQ,nQi, ...). each corresponding to its own 
unique We notice that the dominant C,'s correspond 
to lower harmonics as shown in Table VII. It is interest- 
ing to note that the (0,2,8,34,...) sequence associated 



with the Fibonacci numbers can be considered as the 
second harmonic of the sequence (0, 1, 4, 17, 72, ...) asso- 
ciated with the rational approximants of a^. It should be 
noted that unlike the other two sequences (1, 3, 13, 55...) 
and (5,21,89 ), the sequence consisting of even inte- 
gers (0,2,8,34, ...) does not belong to the sequences as- 
sociated with (7^. However, the fact that it can be re- 
garded as a higher harmonic of one of the sequence 
suggests that the important peaks are determined by 
periodicity and its harmonics, instead of a periodicity. 

Alternatively, the location of the dominant peaks in the 
wave functions can also be related to the number theoret- 
ical properties of cr^ using the Farey tree |Q . Denoting 
the left and right branches of the Farey tree by and 
1, the rational approximants of <t^ can be represented 
by the path 000111100001111... which we will denote as 
O'^l^O^l''.... We conjecture that the dominant peaks in 
the wave functions occur on the sites which correspond 
to the denominators of the rational numbers on this main 
path and on the side paths that differ minimally from the 
main path. 



VI. DECIMATION SCHEME BASED ON 
ARBITRARY FAREY PATHS 

The infinite set of scaling rations can be obtained by 
generalizing the Fibonacci decimation scheme. As de- 
scribed in the previous section, every site in the spin 
chain can be labelled by a symbol sequence correspond- 
ing to a definate path in the Farey tree. We develop a 
decimation scheme based on this. The basic idea under- 
lying the decimation can be understood by considering 
three Farey levels denoted as / (father), m (mother), d 
(daughter). Our decimation equations are, 

«'(« + Qd) = c„(i)*(i + Qm) + ^^„(^)*(^) (25) 
*(« + Qd) - a„(i)^'(z + Qf) + 6„(^)*(^) (26) 

where Qd = Qm + Qf- In order to obtain the recur- 
sion relations for the decimation functions Cn, d„, a„,bn, 
we consider the next level denoted as g (granddaughter). 
Now, Qg = Qf + Qd OT Qg = Qm + Qd depending upon 
whether the Qg is obtained by taking a step in the same 
direction or the different direction from the one from Qm 
to Qd- The above two routes to g level can be denoted 
as 0(1) 0(1) and 0(1) ^ 1(0). 

The recursion relations for 0(1) 0(1) are 

an+i{i) = Cn{i + Qf)an(i) + dn{i + Qf) (27) 

bn+i{i) = Cn{i + Q f)bn{i) (28) 

Cn+i{i) = Cn{i + Qf)+dn{i + Qf)a-^{i) (29) 

dn+i{i) = -dn{i + Q f)bn{i)a;,'^ (i) (30) 



and for 0(1) 1(0) 
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a„+i(i) = a„(i + Qm)c„(i) + 6„(i + (31) 

bn+lii) + Qrn)dn{i) (32) 

Cn+l(i) = a«(i + Qm) + bn{i + Qiyi)Cn^{i) (33) 
rf«+l(«) = -6n(« + <3m)rfn(«)c;;^(«) (34) 

It should be noted that the Fibonacci decimation corre- 
sponds to the symbol sequence 101010... in which case the 
recursion relations can be written in the form in which 
a„ and &„ become redundant. 

In order to explain the infinite set of C's we first need 
the decimation functions along the "main" Farey route 
to which defines four universal scaling ratios (there 
is an 8-cycle in the decimation functions). Other uni- 
versal C's are then obtained by determining the deci- 
mation functions along finite side paths from the main 
route. It should be noted that the most dominant peak 
in the wave functions (corresponding to largest C) occurs 
at 1, 4, 17, 72, ... (see Figs. 2 and 3) and can be obtained 
only by decimation. In the previous RG scheme fl^ ] 
the scaling ratio C corresponding to this sequence was not 
predicted. 

VII. RENORMALIZATION GROUP EQUATION 

The universal functions and the scaling ratios can in 
principle be obtained from RG equations. We have al- 
ready derived recursion relations for the decimation func- 
tions or matrices which can be used as the starting point 
for the definition of a suitable renormalization operator. 
However, in order to get the final equations in a solv- 
able form it is essential to get rid of the discrete lat- 
tice index i. Because the QP potential has been con- 
structed using a period- 1 function, the decimation func- 
tions/matrices an{i),bn{i),Cn{i),dn{i) can be thought of 
being functions of the fractional part of ia, denoted by 
< ia >, only. This gives us a continuous variable. But 
we also need to consider the scaling. Assuming a p-cycle 
for the decimation functions at i = 0, we observe that if 
a„, 6„, c„, dn are defined as functions of the renormalized 
variable x = (— cr)"" < ia >, any decimation function of 
the level n maps roughly onto the corresponding function 
of the level n + p ior x E [0, 1] (see Figs. 5 and 6). This 
indicates that x provides us with a proper continuous 
variable. 

The next step is to rewrite the recursion equations for 
the decimation matrices/functions in a form which uses 
the above renormalized variable x. For simplicity, we re- 
strict ourselves to the Fibonacci decimation. By using 
the relation F„cr = Fn-i — (— cr)", the nonlocal recur- 
sion relations (13-14) can be transformed into local RG 
equations 

Cn+i{x) = c„{-ax - l)c„_i((T^a; + a) 

-d^^(-CTa;)d„+i(x) (35) 



dn+i{x) = -dn{-crx)[dn{~crx -1) + 

Cni-<7X - l)(i„_i(<T^a; + a-)]c^^(-<Ta;) (36) 

This defines our renormalization operator which can be 
now applied to determining the universal cycles. We tried 
that in the case of the critical Harper model (£■ = 0). Our 
detailed numerical studies showed that (c„(a;), c?„(x)) ap- 
proached a six-cycle (ci,d*) {cl,dl) — > (03,^3) 
{-cl,~dl) ^ {c*,-d;) -> i-cldl) -> (cl.dl) where 
all the functions of the six-cycle appeared bounded and 
smooth for the full messure of x values. This suggested 
the functions c* , d* could be expanded in a power series 
and one could use the Newton method to solve for the 
coefficients. However, simple circular complex domains 
centered on the real axis did not give convergent results. 
We believe that the equations can be solved using more 
elaborate domains. It is probably necessary to let the 
centers of at least some of the complex disks to be lo- 
cated off the real axis so that the radius of convergence 
can be increased. This appears true especially for the ex- 
pansion of the inverse of c„ for which this radius is only 
unity if the center of the corresponding domain lies at the 
origin of the complex plain. Therefore, the determination 
of scaling properties using eigenvalues of linearized RG 
equations about the limit cycle remains open. 

As to the RG for an arbitrary Farey path, the main 
difficulty arises from the number-theoretical problem of 
generalizing the relation giving the fractional part of Fna 
to products of the type Q/tr and QmC where Qf and Qm 
are arbitrary integers appearing in the Farey tree. All 
we can prove at the moment is that with the recursion 
(23) the fractional part of Qkcr (with arbitrary Qo,Qi) 
approaches either or 1/2 and the rate of convergence is 
given by (— u)^. However, there are still a lot of details 
which need to be worked out in the general case. 

VIII. SUMMARY OF THE UNIVERSALITY 
CLASSES 

Based on all the infinite set of scaling ratios, we con- 
jecture that the C phase of the model is governed by 
four different fixed points of the renormalization group 
defining four universality classes, which we refer as spec- 
tral, conformal, subconformal and Harper. We show that 
the anisotropy is a relevant parameter as an infinitesimal 
anisotropy takes the system to one of the three possible 
new universality classes, away from the Harper univer- 
sality describing the isotropic limit of the model. Except 
for the conformal universality, the other three universal- 
ity classes exhibit an asymptotic reflection symmetry of 
the wave function. 

We checked in all the above four cases that the uni- 
versality is rigid with respect to replacing the single har- 
monic hi by an odd two-harmonic function. Both the in- 
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finite set of C's as well as the decimation functions were 
found to be universal. 

We also computed the exponent (3 defined by Ostlund 
et al. HJ] which measures the degree of localization 
around the site where the wave function is peaked. 

-/3=^7^ lim[/n[/^]] (37) 

where 

= ^ E ^» (38) 

The p in the above equation is the period of the limit 
cycle. Unlike the scaling ratio C, no close form expression 
for (3 could be found in terms of the decimation functions. 
(3 was computed directly from the wave function. 

The exponent /3 is in the E phase, —1 in the L phase 
and varies between these two values in the C phase. Ta- 
ble VIII shows (3 in four different universality classes. As 
expected the absolute value of (3 is greater at the con- 
formal fixed point compared to the subconformal case. 
Comparing (3 at the onset of localization in the Harper 
and Ising-like cases leads us to conclude that the LRO 
slightly suppresses the degree of localization. However, 
the degree of localization remains constant inside the C 
phase. Since the subconformal C phase existing in fi- 
nite window in A is described by a single fixed point this 
implies that A is an irrelevant parameter in this phase. 

All the other universality classes except the spectral 
one can be described by scalar decimation functions. 
This made the study of the spectral universality a lot 
more difficult. Additional complications came from the 
facts that Emin was non-zero and the phase factor <f> was 
continuously changing. This explains why the wave func- 
tion and the scaling ratios C a-re known with less preci- 
sion in that case than for the the other three universality 
classes. 

It should be noted that the universality classes ob- 
tained here are for the fractal quasi-particle wavefunc- 
tions and do not tell us anything about the scaling prop- 
erties of the fractal quasi-particle energy spectrum. Very 
recently, Chaves |Q showed that the energy level statis- 
tics follow inverse power law distribution with the ex- 
ponent 3/2 throughout the fat C phase. This implies 
that the fractal dimension of the C phase is a constant 
throughout the phase and is equal to .5. p5| . However, 
the multifractal analysis of the spectrum H showed that 
the /(a) curve describing the scaling properties of 
the energy spectrum is different in the four universality 
classes. Therefore, we conjecture that even the multi- 
fractals describing the quasi-particle energy fall within 
the same four classes. We believe that the exact RG 
treatment describing the spectrum as a whole as pro- 
posed by Wilkinson could be generalized to confirm 
this conjecture. 



IX. CONCLUSIONS 

The work described here confirms that the anisotropic 
spin model describing the perturbed Harper with broken 
U{1) symmetry fattens the critical point of the Harper 
model. In analogy with the mode-locked phase, this fat 
C phase existing in a finite window can be described 
as C — locked phase whose scaling properties fall within 
four different universality classes. The conformal and the 
spectral universality classes define the boundaries of the 
phase diagram which sandwich the subconformal class. 
In the isotropic limit, the three universality classes degen- 
erate to the Harper universality. It is rather interesting 
that the subconformal and conformal parts of the phase 
diagram belong to different universality classes. This pro- 
vides a novel example of spectral and magnetic interplay 
showing how a magnetic phase transition (which occurs 
also in a periodic model) affects the critical phase char- 
acteristics of aperiodic systems. The spectral transition 
which exists only in an aperiodic system (and does not af- 
fect the magnetic properties such as the long-range corre- 
lation function) defines its own universality class distinct 
from the conformal and the subconformal cases. 

The E phase of the quasiperiodic models is the KAM 
phase where the wave functions are extended and hence 
the quasiperiodicity has minimal effects. As shown in the 
Fig. 7, in contrast with the C phase, the Hull function 
is smooth here. It is rather surprising that the universal- 
ity in this phase is almost independent of the anisotropy 
g. For the periodic system, the anisotropy is a relevant 
parameter as it provides a preferred axis for spin align- 
ment. The fact that the universal functions of the E 
phase are insensitive to the anisotropy adds a new di- 
mension to the KAM theory: the smooth Hull function 
of the integrable limit A = survives not only the A per- 
turbation in the E phase, but the anisotropy leaves the 
universaliry class unaltered. This is different from the 
C phase where an infinitesimal anisotropy alters the uni- 
versal characteristics. Therefore, we conclude that the 
anisotropy responsible for the long range order does not 
affect the KAM phase and the quasiperiodicity respon- 
sible for the spectral transition does not affect the long 
range order. 

It is interesting to compare our results with those ob- 
tained with a TBM related to the stability of the ground 
state in the classical Frenkel-Kontorova model at the on- 
set of pinning-depinning transition. JTot This pinning- 
depinning transition corresponds to the breakup of KAM 
tori in area preserving maps. In this case the analog 
of the fractal wave function in the C phase is the deriva- 
tive of the Hull function describing a QP critical ground 
state (Fig. 8). The reason why here the wave function 
corresponds to the derivative of the Hull function rather 
than the Hull function itself is the fact that the TBM 
follows from a variation of the ground state. With the 
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initial conditions on the dominant symmetry line (this is 
analogous to our choice of the phase (j)), the Fibonacci 
decimation leads to an asymptotic universal fixed point. 
psf It is interesting to speculate the existence of a C- 
locked critical Hull function in some generalized classical 
Frenkel-Kontorova models or the corresponding area pre- 
serving maps. 

Very recently, Faddeev et al. have shown that 
the Harper equation describing isotropic limit of the XY 
model can be solved using generalized Bethe Ansatz. It 
will be interesting to know if the coupled equations (4-5) 
describing the anisotropic model are amendable to sim- 
ilar treatment. Since in the case of a periodic system 
(cr = 0), the anisotropy introduces an integrable pertur- 
bation on the isotropic model, it is possible that even 
with quasiperiodicity that may be the case. We hope 
that our studies will motivate rigorous analysis of the 
quantum spin models discussed here. 
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FIG. 1. Phase diagram for the quasi-particle excitations: 
With Jx ~ 1 — 5 and Jy — 1 + g, the figure shows the E, 
C and L phases in X — g plane. In the E phase, wave func- 
tions are extended and hence resemble those of a periodic 
system. This phase is also called the KAM phase as the Hull 
functions are smooth in this case. In the C phase the wave 
functions are self-similar (see Figs. 2 and 3) exhibiting atmost 
algebraic localization while in the L phase, they are exponen- 
tially localized. The onsets of E-C and C-L transitions are 
respectively the spectral line (thin line) and the conformal 
line (thick line) . The two lines meet in the isotropic (Harper) 
limit. The g = —1 describes the Ising limit with no E or KAM 
type phase. 
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FIG. 2. (a-e) show the absolute values of the wave func- 
tions, along the line in the phase diagram where the model is 

conformally invariant {Emin = 0). For fixed = 1 > Jy > Q, 
the conformal line is obtained by setting A = 1 and varying 
Jy (see ref. 5): (a) Jj, = , (b) Jy = .99, (c) Jy = .999, (d) 
Jy = .9999, and {e)Jx = Jy = 1 . The figures (a) and (e) 
correspond respectively to the Ising and isotropic limits while 
(b),(c) and (d) arc Ising like. Wo sec that the wave functions 
are highly asymmetric about the peaks in the Ising and Ising 
like cases. However, in the isotropic limit, the wave functions 
recover the symmetry. 



TABLE I. The limiting subcritical (E phase) decimation 
functions for the general anisotropic model {g ^ 0) and the 
Harper limit {g — 0). As the functions do not depend on 
the phase (p, there cannot be any dependence on the lattice 
site i either. For the general anisotropic case, the cycles are 
seen only after diagonalizing the decimation matrices. For the 
Harper limit Emin = —Emax while for the anisotropic model 
Emin is positive. Note that Emin of the anisotropic case goes 
into i5 = of the Harper model as g goes to zero. 



E c* d" cycle 



FIG. 3. (a) shows the wave function below the conformally 
invariant point for the Ising model (Jy = 0) with A = .95. (b) 

is the wave function of the critical model with no diagonal 
disorder, (c) shows the wave function at the onset of the 



spectral transition with 



1, Jv 



1.5, and A 



1. In 



all these cases, the wave function is always asymptotically 
symmetrical (about i = 0). In the limit A — » 0, but finite, the 
wave function is fully symmetric (b). 



FIG. 4. (a-d) respectively show the Hull function |'!/'(a;)| 
vs X =< ia > in the four universality classes, namely 
the Harper [E = 0), conformal [Jy = 0, = 1, 

\ = 1, E = Q), subconformal (critical model with A = .5, 
Emin = -1.420111920492815), and spectral (J^ = 1, 
Jj, = 1.5, A = 1, Emin = 0.17429863452742). 



9^0- 

F, 

Emax 

5 = 0: 




1, 0, oo 



0, 1, oo 



1, 0, — oo, —1, 0, — oo 



a-\-a-\-a-^ 



0,-1,-00,0,1,00 6 
u, —a, a 3 
-a 1 



FIG. 5. (a-d) show the decimation function Cn{x) 
(n = 6 — 9) for the critical Harper model with E = 0,4> = 1/4. 
As a consequence of the 3-cyclo for the absolute values of the 
decimation function, ce{x) = ±09(0;). 



TABLE II. Estimates for Emax of the Harper model with 
A = .5 obtained from the condition Cn(0) = a~^. The n here 
refers to the n*^ order Fibonacci number with F2B = 46368. 



n E 



FIG. 6. (a-d) show the decimation function Cn(x) 
(n = 6 — 9) for the Ising model with E = 0,(j) = 1/4 at 
the conformally invariant point. As a consequence of the 
3-cycle for the absolute values of the decimation function, 
ceix) = ±C9{x). 



FIG. 7. The wave function (a) and the smooth Hull func- 
tion (b) in the E phase at A = .5, Jx = 1, Jy = 1.5, and 
E = Emax = 2.620638044392. UnUke the C phase, the Hull 
function is smooth. 



10 


2.144122230954376 


13 


2.144103497103577 


16 


2.144103742043903 


19 


2.144103738794523 


22 


2.144103738837256 


25 


2.144103738836694 



FIG. 8. The derivative of the Hull function at the onset of 
breaking of analyticity in the Frenkel-Kontorova model. 
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TABLE III. Ejnin for the Ising model obtained by the 
diagonalization (size N finite) and the decimation scheme 
(AT = oo). 



A N 



0.95 


1597 


1.51567E-2 




2584 


1.51484E-2 




4181 


1.51481E-2 




oo 


1.514806760988E-2 ± l.E-14 


0.995 


4181 


6.1929E-4 




6765 


6.1862E-4 




10946 


6.1853E-4 




oo 


6.18420285E-4 ± l.E-12 



TABLE IV. The decimation functions at i = for the 
Ising model with A = l, i? = 0, ^ = 1/4. The table clearly 
shows the asymptotic six-cycle. 



n c„(0) d„(0) 



3 


-2.245709795011781 


-0.6002603574022150 


4 


8.109008292790231 


-1.272168729380112 


5 


1.010919961015582 


-0.3666408658531275 


6 


3.867229946686237 


-0.7949593051332248 


7 


-7.553298439810634 


1.182294059162462 


8 


1.100178678647771 


0.3764230488289012 


9 


-3.710289833700136 


-0.7701746581532440 


10 


7.758396265078099 


-1.194676151297624 


11 


1.085391362405107 


-0.3738057364522718 


12 


3.746636404387838 


-0.7734760239494146 


13 


-7.718184336265150 


1.192015324586235 


14 


1.089002996132289 


0.3742871587654810 


15 


-3.738988672892117 


-0.7728088313940935 


16 


7.727825304953185 


-1.192588926144134 


17 


1.088173364711348 


-0.3741752202327470 


18 


3.740771704903936 


-0.7729526819664260 


19 


-7.725577268094126 


1.192454968911469 


20 


1.088369359725458 


0.3742010466453257 


21 


-3.740355675780264 


-0.7729194653703071 


22 


7.726105073162105 


-1.192485868259755 


23 


1.088319511696401 


-0.3741951424798995 


24 


3.740445328629009 


-0.7729293725799985 


25 


-7.726061715596284 


1.192499566226364 



TABLE V. Estimates for the conformally invariant point 
of the Ising model where the third harmonic with a = .2 has 
been added to h{i). The estimates were obtained by impos- 
ing the condition c„(0) = 1.08832 for </> = 1/4 (see the first 
column corresponding to n = 23 in Table IV) . 



n A 



11 


1.367474785 


14 


1.366202286 


17 


1.366282249 


20 


1.366277923 


23 


1.366278173 


26 


1.366278174 



TABLE VI. ("'s corresponding to some of the dominant 
peaks in the Harper, Ising (both conformal and subconfor- 
mal), and the spectral universality classes. 



sequence Harper conformal subconformal spectral 



0,1,4,17,72,... 


.516 


.744 


.676 


.72 


1,1,5,21,89,... 


.315 


.622 


.590 


.54 


1,2,9,38,... 


.256 


.297 


.337 


.58 


1,3,13,55,... 


.143 


.496 


.231 


.37 



TABLE VII. C's corresponding to a sequence and its har- 
monics for the Harper model. 



harmonic sequence C, 



1 


3, 13, 55, 233,... 


0.143 


2 


6, 26, 110, 466,... 


0.1456 


3 


9, 39, 165, 699,... 


0.1485 


4 


12, 52, 220, 932,... 


0.0183 


5 


15, 65, 275, 1165,... 


0.025 
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TABLE VIII. The exponent /3 for the four universahty 
classes. 



Universahty Class 



Harper -.639 ± .005 

conformal -.602 ± .001 

subconformal -.423 ± .001 

spectral -.175 ± .005 
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